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The Skornyakov-Ter-Martirosyan three-boson integral equations in momentum space are transformed into 
differential equations. This allows us to take into account quite directly the Danilov condition providing self- 
adjointness of the underlying three-body Hamiltonian with zero-range pair interactions. For the helium trimer 
the numerical solutions of the resulting differential equations are compared with those of the Faddeev-type AGS 
equations. 
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I. INTRODUCTION 

The discovery of a weekly bound state of two helium 
atoms Q| and problems concerning the stability of Bose con- 
densates of alkali atoms (see, e.g., |2[) stimulated many calcu- 
lations of the properties of three-particle systems determined 
by pair interactions with a large scattering length oq compared 
to the range ro of the pair forces, 



£0 



< 1. 



(1) 



On the one hand, this condition poses problems in the nu- 
merical solution of the Faddeev equations with realistic atom- 
atom potentials, on the other hand it forms the basis of the 
zero-range model of the two-body interaction ((ro — > 0)). This 
model determines the motion of particles beyond the pair- 
interaction region, admits an analytic simplification, and can 
be used for the description of real physical systems provided 
the condition Q is fulfilled. 

The zero-range model (ZRM) for three-body systems, 
though having obvious advantages and wide application, has a 
considerable drawback. The Hamiltonian of these systems is 
not self-adjoint (see, for example, |3] ) and the Schrodinger 
equation has quadratically integrable solutions at any en- 
ergy. This fact was pointed out in |4] in analyzing poor at- 
tempts to solve the nd-scattering problem with the use of the 
Skornyakov-Ter-Martirosyan (STM) integral equations |5] 
within the zero-range model of nucleon-nucleon interaction. 
In particular, at large momenta k of the relative motion of a 
particle and a pair the asymptotic behavior of the wave func- 
tion was shown to have at any energy Z the form 



,, 7 , . sinQornlT)) cos(/j ln(fc)) , 
y{k,Z)=A - 2 + B r? +o( 1J ). (2) 
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Here, the constant /jo depends on the ratio of the particle 
masses, and only the coefficients A and B depend on the en- 
ergy. To determine the three-body spectrum, Danilov 0] sug- 
gested to use the relation between the coefficients A(Z) and 
B{Z) with an energy-independent parameter y, 



A(Z)=yB(Z), 



(3) 
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following from the orthogonality condition of eigenfunctions. 
In the same year Minlos and Faddeev |3] showed that the 
Danilov condition is a special case among possible extensions 
of the Hamiltonian to a self-adjoint one. Even after the ex- 
tension, a solution to the STM equation has one free param- 
eter. This opens the possibility of describing real three-body 
systems by adjusting the free parameter to a known spectral 
point 0|. 

Unfortunately, the Danilov condition for an unambiguous 
solution of the STM equations practically cannot be used in 
numerical calculations. But it was used for an analytic inves- 
tigation of the ZRM three-body spectrum 1 6 ] . In this reference 
it has been shown that a three-body collapse (the Thomas ef- 
fect 1 7]) is a specific feature of the ZRM. 

There are several approaches using the ZRM beyond the 
scope of the STM equations, among them the adiabatic ex- 
pansion in configuration space (see, e.g., 0]). Just in this 
approach the so-called Efimov effect was observed, i.e., the 
fact that three-body spectra concentrate on zero total energy if 
the two-body scattering length |ao| tends to infinity |9|]. The 
result of adiabatic expansions is an infinite system of differ- 
ential equations coupled in terms of first derivatives. The 
problem of non-self-adjointness of the three-body Hamilto- 
nian is solved by cutting off the effective interaction at small 
distances. This method of regularization implicitly introduces 
three-body forces. In this case, the cut-off radius plays the 
role of the free parameter. 

Three-body forces are introduced more explicitly in the ef- 
fective field theory (EFT) II 1(11 . causing an integral equation 
that is similar to the STM equation, but contains artificial 
terms. A free parameter enters the phenomenological terms 



2 



of the kernel of the integral equation and the free terms of the 
integral equation for the scattering problem 1 10]. 

Finally, we would like to draw attention to the two-pole t- 
matrix model for the description of three-boson systems lUlll 
which employs the Alt-Grassberger-Sandhas (AGS) 1 12J ver- 
sion of the Faddeev equations [13]. The position of the sec- 
ond pole on the unphysical sheet is treated in this model as a 
parameter of the integral equation. When the position of the 
second pole tends to infinity, the STM equations are repro- 
duced 1 1 1 ]. This method deals with compact equations and is 
well suited for numerical calculations, but is of little use for 
analytic considerations. 



II. FORMALISM 

In the present paper we follow another strategy. Transform- 
ing the STM equations into differential equations allows us to 
take into account the Danilov condition quite directly. For this 
purpose, we construct an infinite system of differential equa- 
tions in momentum space with a very simple relation between 
the equations. 

The homogeneous part of the STM equation for the elastic 
scattering amplitude |5] of a boson of mass m on a two-boson 
bound state of energy £ = —x 2 /m can be represented in a form 
convenient for further analysis, 

J®* (4) 



nj k 2 +k 2 -k t k + X 2 ^/X 2 + k 2 3/4-. 

where X 2 = —mZ. The function F(k) is related to the 
wave function of the three-boson system cp(fc) via F{k) = 
k(y/ -mZ + k 2 3/4- x)(p(fc). It should be noted that Eq. (g) 
implies F(—k) = —F(k). 

For further transformations of 10} it is convenient, follow- 
ing H, to substitute the variables k = (F(k{t)) = F(t)) 
and use the Mellin transformation 5F(s) = J °° t s ~ l F(t)dt 

=L(s)?(s) +L(s)2^ J dt j£ '— Fit), (5) 
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where 



8 sin(^) 



y/3 s cos(fs)' 
By inverting (1 — L(s)), Eq.© goes over into 



J( J )=J„( J ) + T ift T 2? (dt ' 



l-L(s) X 
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Here, the general solution of 

(i-z,0?))3^0O = o 



F(t). (6) 



(7) 



for s = ±//jo (/JO = 1 .00623 . . .) had to be added. We recall that 
Eq.Q was used by Minlos and Faddeev |6] to study the STM 
equation spectra. In the variable t the solution is 

F as (t) oc sin(/j liu). 



When carrying out the backward Mellin transformation we 
see that the function F can be written as a sum 



■■FO- 



LK 



(8) 



over the residues of the function (1 — L(s)) . The term Fq 
is generated by the poles at the points s = ±//Jo and the terms 
Fj correspond to the real positive solutions of the equation 1 — 
L{X, •) = 0. It can be shown that the components of the function 
F satisfy the differential equations 



, d d . j 2uq 2y 



D(fJo) ^±l-2f 



F(t) = 0, (9) 



r d d , , 2Xj 2f 

^- Fi) -x 2 Fr 



D{X t ) £±i-2f 



F(t)=0, (10) 



where D(X,-) = (1 -L(s))'\s=x, and iD( m ) = (1 -L(s))'\ s= , m . 
Boundary conditions at t = 1 , (k = 0) follow from the above- 
mentioned antisymmetry of F(k), i.e., Fj(t)\ t= \ =0; i = 
0, 1, .... At large values of the argument one can use either the 
property of boundedness, or the asymptotic form following 
from condition (0: F(?)| t _>oo — * sin(/Jolnf + 8). Comparing 
the Danilov condition Q with this asymptotic form of F we 
get an extra condition for the spectrum 



8 + jUoln(\/3 /X) = const{Z) +nn. 



(11) 



It can be verified that the asymptotics of F coincides with 
that of Fq and the contribution of the remaining terms de- 
creases at infinity. At x/X = the system of equations J91 1 0I > 
becomes uncoupled and has a simple solution, F = Fq = 
sin (/jo In (f )), which gives, together with condition 1^2, the 
spectrum E„ = £bexp(— 2nn). This spectrum contains the 
tree-body collapse at large X for n < |6], and Efimov's con- 
centration towards the point Z = at extremely small x for 
n > 0. At small finite x we can restrict ourselves to Eq. l|9) 
that belongs to a well known class of differential equations 
(Heun) (3. 



III. RESULTS 

Let us first demonstrate the convergence of the solutions for 
an increasing number of equations fl9HUt . To this end, we 
denote the solution of this system of equations by Fr which, 
upon substitution into the right-hand side of ©, gives the 
function Fl- The degree of proximity of Fr and Fi shows to 
what extent the solution of the system fl9H0i is close to that 
of the integral equation. Figure [2 shows this convergence for 
N = 1,2,3. The energy parameter Z/e = (X/x) 2 was chosen 
to be 1.57. This value corresponds to the calculations lf]~5ll 
for the binding energy of a helium trimer. In this way we 
fix the free parameter characteristically occurring in all STM 
treatments. The good convergence achieved already for N = 3 
indicates the efficiency of our differential equations approach. 

Instead of making a similar consideration for the bound- 
state spectrum, we compare our present results with alterna- 
tive calculations. In this context we use the binding energies 
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FIG. 1 : Convergence of the solutions of Eqs. J91 101 for increas- 
ing N: Curve 1 shows Fi, curve 2 shows Fr 
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FIG. 2: Trimer binding energies: 1 calculations by Eqs. J91 101 : 
2 calculations II ill : 3, 4 calculations for three- and four-level 
trimers; 5 calculations 1 15, 16] for a helium trimer 1 15, 16] 



of a helium trimer, obtained in for realistic pair poten- 

tials, and the calculations of the binding energies via Faddeev- 
type AGS integral equations in the framework of the above- 
mentioned two-pole pair f -matrix II ill . Since, depending on 
the position of the f-matrix pole on the unphysical sheet, one 
can obtain as many bound trimer states as one wishes, we will 
label the highest bound state by "1" and the following one by 
"2". 



Figure [2] shows energy E% as a function of energy E\ cal- 
culated with Eqs. A9I1 0I > for N — 21. Also shown are the 
calculations for a two-level trimer within the two-pole pair t- 
matrix [11], and our calculations in the same model when the 
parameters of the pair f-matrix admit the existence of three 
and four bound states of a boson trimer. All the energies are 
given in units of the dimer energy. We see that there is only 
a rather small difference between the corresponding curves. 
Quite interesting is the fact that the calculations for a helium 
trimer II 5l Hal with realistic pair potentials lie just between 
these curves. 



Thus, we have demonstrated that Eqs. J9ll0i represent a 
very efficient tool for calculating three-boson trimers. These 
equations can also be extended to scattering problems. 
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